Theory of Hysteresis Loop in Ferromagnets. 
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We consider three mechanisms of hysteresis phenomena in alternating magnetic field: the domain 
wall motion in a random medium, the nucleation and the retardation of magnetization due to slow 
(critical) fluctuations. We construct quantitative theory for all these processes. The hysteresis is 
characterized by two dynamic threshold fields, by coercive field and by the so-called reversal field. 
Their ratios to the static threshold field is shown to be function of two dimensionless variables 
constituted from the frequency and amplitude of the ac field as well as from some characteristics of 
the magnet. The area and the shape of the hysteresis loop are found. We consider different limiting 
cases in which power dependencies are valid. Numerical simulations show the domain wall formation 
and propagation and confirm the main theoretical predictions. Theory is compared with available 
experimental data. 
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I. INTRODUCTION 

The hysteresis loop (HL) has been first studied more 
than century ago Q. However, the understanding of this 
process in thin magnetic films as well as in bulk magnets 
is still rather poor. Many efforts have been devoted re- 
cently to prediction (see |^-^) and experimental verifica- 
tion (see 1^,^ , ) of the scaling behavior of the hysteresis 
loop area (HLA) as a function of the applied magnetic 
field frequency and amplitude for thin magnetic films (for 
a brief review of HLA scaling results see [^). The scal- 
ing behavior of the HLA have been first reported in the 
pioneer work for 3D magnets. While there exists an 
extended literature on the hysteresis of 3d magnets, the 
properties of HL in 2d systems are much less known. 
There are only few articles devoted to the HL in ultra- 
thin ferromagnetic films [|o|, 0, §], though the 
hysteresis effects have been found as a side effect in many 
others (see, for example Critical exponents 
found in the experiments with thin films vary dramat- 
ically for different materials (see e.g. and ^) and 
probably for different regimes. Different authors disagree 
each to other (see the already cited articles [@, B| and ||]) 
and also disagree with numerical simulations 1 2 ] . 

Several years ago mean field type models with single 
pof or many jl^] relaxation times have been applied to 
analyze experimental data. The authors of as- 
sumed that the HL was controlled by the nucleation pro- 
cess. These authors predicted the logarithmic depen- 
dence of the coercive field he on the rate of the applied 
magnetic field h. In a recent experiment |^ it was found 
that the HLA depends on frequency of applied field as a 



power with a small exponent (^ 0.03 — 0.06) or, possibly, 
logarithmically. However, in the framework of the same 
approach the HLA must behave also as logarithm of h 
with the same coefficient (exponent). This dependence 
has never been observed in the experiment. Therefore, 
we propose a new analysis of such a HL in this article. 

The purpose of this article is to formulate a rather 
general approach to magnetization reversal mechanisms 
and to indicate several important measurable character- 
istics of the HL besides of the HLA. We will see that 
these characteristics are governed by two dimensionless 
parameters combined from the field frequency w, its am- 
plitude and characteristics of the magnetic material. 
Everywhere in what follows we assume that the external 
field varies harmonically in time h{t) = ho sinut. 

The hysteresis behavior may have various origin. It can 
be mediated by nucleation process, by the DW propaga- 
tion or simply by retardation of the magnetization due 
to fluctuations. We consider all these mechanisms and 
establish conditions at which one of them is dominant. 

Defects play an important role in the DW propagation. 
They create a finite threshold value hp of static magnetic 
field necessary for the DW depinning. The threshold 
field hti in the dynamical problem can differ substan- 
tially from hp. We find that in a medium with defects 
the moving DW, passing rare extended defects, may form 
bubbles of reversed spins. These bubbles play an impor- 
tant role as prepared nuclei in the next half-cycle of the 
magnetization reversal. 

In this article we consider magnets of the Ising (uni- 
axial) symmetry. Their properties may be very different 
depending on the strength of the anisotropy. In the ex- 
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perimentally studied films the anisotropy was very weak. 
In this case the domain wall width is large in comparison 
to the lattice constant. On the contrary, in the origi- 
nal Ising model the anisotropy is assumed to be large 
and DW width I is simply the lattice constant. However, 
these different models become equivalent after a simple 
rescaling: the DW width should be accepted as a new 
elementary (cut-off) length. It means that we consider a 
spin cluster of the linear size I as a new elementary spin. 
This approach allows us to apply the Ising model supplied 
with the Glauber dynamics for numerical simulations. 

Peculiarities of the two-dimensional situation are: 
much higher mobility of the DW as well as much stronger 
fluctuations. This makes the experimental situations as 
well as the theoretical description much more diverse 
than those for a 3D magnet. 

This article is organized as follows. In sections ^ and 
III we consider the individual DW motion. Equation of 
motion is formulated and justified in Sect.^ It is solved 
in Sect. [11. In the same section we introduce character- 
istic fields hti, ht2, he, hr, the HLA A, find the scaling 
arguments and analyze several limiting regimes in which 
simple power scaling is valid. In Sect. IV the process of 
the bubble formation is studied. Sect .[V| is devoted to the 
HL controlled by the nucleation process. The HL driven 
by strong magnetization fluctuati ons , especially near the 
Curie point, is considered in Sect, VI, Numerical simula- 



tions of the HL and the domain structures for the 2d Ising 
model with Glauber dynamics which supports results ob- 
tainedin Sect. Ill, Sect. IV are presented in Sect. VII, In 



Sect. VIII we summarize our results and compare them to 
the experimental data. In the rest of this article we use 
the notation h for magnetic field, m for magnetization 
and M for the total magnetic moment of a magnet. 



II. DOMAIN WALL MOTION IN A RANDOM 
MEDIUM 

As we already mentioned in the Introduction, our 
starting point is an impure ferromagnet with either weak 
or strong Ising anisotropy. The soft spin version of the 
system is then given by a model with a bare domain 
wall width I a, where a denotes the original lattice 
spacing. The imperfections in the model may be in prin- 
ciple either of random bond (i.e. random Tc) or random 
field type. We will argue below, that in the region we are 
mainly interested in, namely above the depinning thresh- 
old, both types of impurities act essentially as random 
field impurities. 

As it was shown by Bausch et al. jl^ (see also ) , 
equation of motion for a domain wall without overhangs 
can be written in a following way 



1 dZ 



rv- +/i + 77(x,z) (1) 



where Z{x,t) denotes the interface position and g = 
1 + (VZ)^. 7 and r are the domain wall mobility and 
stiffness, respectively, h = ^bHM, where H is the ex- 
ternal magnetic field and M is the magnetization. Finally 
77 denotes the random force generated by the impurities. 

For broad domain walls F « J/{a^^^l), where D de- 
notes the dimensionality of the wall. For narrow walls F 
depends in general on J, T, and the disorder strength in 
a complicated way [|l8| . 

The random fields ri{f ~ (x, Z)) generated by imper- 
fections is assumed to be Gaussian distributed and short 
range correlated with r]{'r) = and 



r]{r)r]{r') ^ r]n^+'^6i{r ~ f') 



(2) 



Here Si{f) denotes a delta- function smeared out over a 
distance I. As was first argued by Narayan and Fisher 
in the region above the depinning threshold ran- 
dom bond and random field impurities act in a similar 
way. This can be seen most easily from an example of 
two ratcheted potential (see figure 1), one for random 
bond and another for random field impurities. Although 
the potential V{Z) = ri(x, Z') dZ' in the random field 
case (figure lb) has fluctuations which scale like Z^/^ , it 
leads to the same random forces as the random field po- 
tential (figure la) which shows of order one fluctuations. 
Therefore, we restrict ourselves throughout this paper to 
the case of random field impurities. Moreover, we will 
assume that the disorder is weak, i.e. that the condition 



F>77; 



(3) 



is fulfilled. 

For |V^p^l,.gssl and the equation of motion takes 
the form considered previously jl^ ||l^ ||l9[| fl^. Below 
we summarize some of the results found in [|l9| |14[] : Since 
the disorder is weak (see (||)), the interface is essentially 
fiat on length scales L <^ Lc where 
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rjl 



2/(4-15) 



(4) 



is the so-called Larkin length. On larger scales the wall 
can adapt to the disorder and, as a result, it gets pinned 
for driving fields h ^ hp with 



D/{4-D) 



< 77 



(5) 



for the pinning threshold. If h exceeds hp, the wall starts 
to move. For h^ hp the influence of the disorder is weak 
and the velocity is proportional to the driving field 



D 



7/1 



(6) 



Corrections to this relation can be considered in the 
framework of high-velocity expansion which can be ex- 
pressed as a power series in 
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Lc 



hp-f 



1/(^-0 



(7) 



where is dynamical correlation length which diverges 
as w — !■ 0; z and C are the dynamical and the roughness 
exponent, respectively. Outside the dynamical critical 
region, i.e., for h — hp ^ hp, z = 2 and ( = 0, respec- 
tively. 

If the driving field h is so small, that ^t, ^ Lc, the 
high-velocity expansion breaks down and a (functional) 
renormalization group calculation has to be applied [p^ . 
This leads to a renormalization of the mobility constant 
7 7cff with 

(l/3)(4-D-C) 



7cff 



(8) 



After integrating out the interface fluctuations on the 
length scales L ~ the effective equation of motion 
for the interface profile Z{x,t) — (Z(x, t))^ ^ on large 
scales is given by 

1 dZ 

-J- ^rV^Z + h-hp + T],s{^,vt) . (9) 

Here ()j ^ denotes the spatial and time average over 
scales S,y and ty, respectively and rjcs is the renormalized 
random field which acts as a thermal noise. Since the lat- 
ter leads to an interface roughness characterized by the 
exponent = (2 — D)/2, we may neglect the influence 
of the random field on these length scales. The mean 
velocity of the interface is given by 



h — hr 



h-h 



^«1 



(10) 



where 9 — {z — C,)/{2 ~ C,). z and C take now non- 
trivial values, which can be calculated in by e-expansion 
in D = 4 — e or determined numerically. For D = \ the 
e-expansion gives C = 1 and z = 4/3 and hence = 1/3 
whereas the numerical values are C, — 5/4, z k, 1.42 and 
« 1/4 To unify our results for the domain wall 

velocity we rewrite (^), ( p^ ) as 



X , a; 3> 1 



(11) 



On the length scales Lc ^ L <^ ^y the domain wall is 
rough 

w{L)^({Z{K,)~Z{^,)fy/' ^ 1{L I Lcf . {12) 

\ I |Xi-X2|=l 

Bumps in the domain walls which emerge from random 
clusters heal on time scales 



I fh 



V \v\ 



-fhp \hp 



h 

-. 1 



(13) 



However, on larger scales L £^y the random field acts 
merely as a thermal noise and the roughness exponent is 
reduced to C = 1/2 and C = (log) in D = 1 and £) = 2 
dimensions, respectively. 



III. MOTION OF A RECTILINEAR DOMAIN 
WALL. 



We start with description of a rectilinear domain wall 
motion. In previous section we demonstrated that the 
domain wall roughness can be ignored on a time scale 
t > ty. Thus, locally the domain wall moves as a straight 
line. In some experiments only one domain wall survives 
(see sect. 7). In this case the model problem of rectilinear 
domain wall motion is close to reality. In other cases this 
problem is an important part of more complex problem 
describing either local properties of domain wall motion 
or the order of magnitudes for the global motion. Thus, 
we consider the motion of a rectilinear domain wall under 
the action of magnetic field antiparallel to the magneti- 
zation. 

We have mentioned already that the fluctuation bend- 
ing of the domain wall can be neglected if the charac- 
teristic time of the process is much more than the bump 
healing time ty cx . This requirement suggests 

that ujty <^ 1. Anyway, this requirement must be satis- 
fied since, otherwise the average position of the domain 
wall almost does not change during half a period of os- 
cillations. 

The domain wall is assumed to be fixed at the left 
boundary of the sample Z = at the initial moment. 
We will solve equation of motion (^ for domain wall 
coordinate for harmonically oscillating magnetic field 
h — hosaiLot. Instead of integrating it over time, we 
integrate it over field by a following change of coordi- 
nates: 



dt 



dh 



(14) 



After integration we find an expression for Z vs. mag- 
netic field h: 



Z 



I 



h — hr 



^ Jh, \ hp J - /i2 



dh 



(15) 



This equation is correct for h > hp. For smaller value of 
h the domain wall docs not move: Z — const. Equation 
jlSl ) should be complemented by a prescription to change 
the sign of the square root each time as h reaches its max- 
imum or minimum value ±ft,o and by an initial condition 
Z = Oatt — h — 0. The second necessary prescription is 
to substitute h — hp hy —h — hp when h is negative. To 
transfer from the coordinate Z to the magnetic moment 
M , a following representation is useful: 
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M = M. 



2Z-L 



(16) 



where is the saturation magnetic moment. First of all 
we find two important boundary values for the amplitude 
ho which separate hysteresis loops of different shapes. 
The first of them is the dynamic threshold field hti, the 
smallest value of ho at which the domain walls reaches 
the right boundary of the sample Z = L. At hg < hti 
the magnetization is not reversed fully and the hystere- 
sis loop is asymmetric (see Fig ||a). At larger values 
of ho the hysteresis loop is symmetric under inversion 
h — > ~h,M — > -M (Fig. |c). The value of hn is 
determined by a following equation: 



ujL 
2jhp 



hti /hp 



dx 



(17) 



Note, that at ho = hti the hysteresis loop is symmet- 
ric with respect to reflections in the axis h and A4 (Fig. 
|b). From Eqn. (|l^ it is seen that the ratio hti /hp 
is the function of one dimensionless variable ujL/('jhp). 
This scaling relationship will be analyzed in details later. 

The next notorious field is ht2 defined as the value of 
ho at which the domain wall reaches the right end of the 
sample Z = L during one fourth of period, just at h — ho- 
Equation which defines ht2 is rather similar to that for 
hti- 



ht2/hp 



/(a^-l) 



dx 



Viht2/hp)^ - x2 



(18) 



It differs from Eq. ( p^ by the absence of factor 2 in 
denominator of the l.-h. s. The hysteresis loops corre- 
sponding to ho > ht2 acquire characteristic "whiskers" , 
single-valued pieces of the curve M{h) which are absent 
in hysteresis curves for ho < ht2 (see Fig. ||d). 

At a fixed ho > hti it is possible to define the coer- 
cive field he by the requirement M{hc) = 0. Finally, for 
ho > hti the so-called reversal field hr can be defined 
as a value of the field h at which the magnetic moment 
reverses fully. At the value h — hr two branches of the 
hysteresis curve intersect each other. In other words, 
Z{hr) — L. At h between hr and ho the magnetic mo- 
ment remains a constant M — M^. The values he and hr 
are shown in Fig. ^. Using equations (|l^) and (|l^), we 
find: 



LOL 



Equation for hr reads 



J ( h-hp 



dh 



(19) 



ljL 
ihp 

ihp 



^ (h-hp 



dh 



{ho>ht2) (20) 



(21) 



h — hr 



dh 



h^ 



{hti <ho< ht2). (22) 



The ratios hc/hp and hr/hp are functions of 2 dimen- 
sionless variables u = ujL/["fhp) and v = ho /hp. Note 
that by knowledge of W2 — ht2/hp (or wi = hti/ hp) as 
a function of parameter u = ujL/jhp one can restore the 
function f{x) solving the Abelian equation: 



2 d 

TT dx 



u{w)wdw 

TP 



(23) 



The area of hysteresis loop A can be also expressed in 
the integral form: 



Ah 

A ^ ^hpMs - 4— 
L 



Z{h)dh. 



(24) 



Now we proceed to the analysis of the hysteresis loop 
characteristics. It is controlled by parameters u = 
bjL/^hp and v — ho/ hp. We start with small u ^ 1 
imposing no restrictions on the value v. First we show 
that, at small u, the fields hti and ht2 are close to hp. 
Indeed, it is clearly seen from eqns. (|l|) and (|l|). Solv- 
ing them approximately and employing the asymptotic 
formula for f{x) at small x (pTl), we find 



1 rhti/hp 

'"-Til <^-" 



dx 



^J{hti/hp) - X 



(25) 



x-l 



Introducing a new integration variable s — (Ti^yTTTT^' 
eqn. (|2^) can be transformed as follows: 



u 2 = —;=[-, 1 

V2\hp 



e+i/2 



5(0 + 1,1/2), (26) 



where i?(x, y) is the Euler beta-function. Using its stan- 
dard representation pit], one finds: 



hti 
hr, 



(l + l/26l)r(6'+l/2) 



27rr(6') 



(27) 



In a similar way the asymptotic of the field ht2 can be 
established. The ratio {ht2 — hp) /{hti 
depend on the parameter u if u is small 



hp) does not 



ht2 - hp 
hti — hp 



— 22» + i . 



(28) 



The shape of the hysteresis loop and its area depends 
not only on the parameter u, but also on the parameter 
V = ho/hp. If V is close to 1, i.e. if ho is also close to hp, 
one can employ the low-field asymptotic of the function 
f{x) ( pT| ) also in equation ( pO| ) for hr- Then, for ho > ht2 
approximate equation for hr reads: 
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1 



hr/hp 



dx 



(29) 



Introducing a new integration variable s ~ , we find 



{v-lf+^/^B{e + l,l/2;w); 



hr 



ho ^ hp' 



(30) 



where B{x, y; w) is incomplete _B-function, defined by an 
integral: 



/•W 

B{x,y;w) = / 3^-^(1 - s)y-^ds. 
Jo 



(31) 



/o 

hr — h 



ho — hp 



- in this limit is a function of 



Thus, the ratio w 

only one variable u/{v — 1)^+^/^. For Hq close to ht2, the 
field hr is also close to hp. If u/{v — 1)^+1/^ ^ the 
ratio w becomes small: 



V2v 



'+l)(v-l)''+l/2 



(32) 



For completeness we present here equation for hr in the 
range hti < ho < ht2 without derivation: 

u ^ {v- l)»+i/2 [25(61 + 1, 1/2) - B{9 + 1,1/2; w)] . 

(33) 

If u <C 1, but w ^ 1, the value hr depends on the product 
uv. If it is small, then hr is still close to hr, 



ip. 



[(0 + i)H'/^'+" 



l)w/ioL 



(34) 



In the opposite limiting case {u ^ l,uw ^ 1) the field 
hr is much larger than /ip, and the asymptotic f{x) ~ x^ 
is not valid. In this case equation for hr reads: 



/ ujL ( , ujL 

2/io 

7 V 7 



(35) 



Still the reversal field hr is much less than the amplitude 
/iQ- For coercive field /ic in the range u <C 1, w 3> 1 one 
finds: 

^-1 = f^-l)-2-V(m). (36) 
/ip \hp J 



he 

hr, 



1 

71' 



(37) 



The shape and the area of the hysteresis loop depends 
on the same parameters u — LjL/(^hp) and v — ho/hp. 
We consider the range of small u ^ 1. If w is close to 
one, the low-field approximation for function f{x) may 



be used in equation (|lq)for z. Together with equation 
(|^) it results in a following relation: 



M = M, 



B [0 + 1,1/2,^^ 



B(d + l,l/2, 



hr — hr 



(38) 



For — hr < h < hp on the lower branch of the hystere- 
sis curve and for ^hp < h < hr on the upper branch 
M = zLMs- We presented here only a comparatively 
simple case /iq > ^t2- The case /iq < ^t2 is more compli- 
cated and we do not present here the explicit formulae 
for magnetization in this range of amplitudes. 

For u <^ 1, V ^ 1 and m?; <C 1 we find in a similar way: 



2 1 h-hp 



(hp < \h\ < hr). 



(39) 

and M = Ms for /i^ < |^| < /lo- Note that the hys- 
teresis loop in this range of parameters is narrow, i.e. 
hr, he <C /lo- 
in the interval u <C 1, uw 3> 1 the value of hr, accord- 
ing to eqn (p5|), is much larger than hp. One can neglect 
the threshold field and employ the linear asymptotic for 
f{x) solving eqn. (15). The result is: 



Z{h) = ^iho- \/hl-h^). 



and 



M = M, 



1 - 2^(/io - Jhl-h^) 



(40) 



(41) 



The latter two equations are valid in the interval of 
fields \h\ < hr. Beyond of this interval the total mag- 
netic moment is fixed at its saturation value. Note 
that a new characteristic length appears in the problem 
= ^hQ/uj. Below we present results for the hysteresis 
loop area A in the range of small u for different values of 
V without repeating of analogous calculations: 
i) w - K 1: 



A « AMshp 



1 + ^-^ '- / B{e + l,l/2;x)dx 



ii) w S> l,uv <^ 1 : 



A « '^Mshp 



e + 2 



(42) 
(43) 

(44) 



iii) 3> 1: 

In this case hr hp. Therefore all essential results for- 
mally coincide with those for u 3> 1 and the reader is 
referred to a corresponding subsection. Note that in the 
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cases i) and ii) A is close to a constant value = ^Mshp, 
but the deviation A — Aq scales with two parameters u 
and V. 

Now we proceed to a simpler case of large u 3> 1. In 
this case the dynamic threshold field hti is much larger 



than hp'. 



hti = 



(45) 



Therefore, one can neglect hp when integrating eqn ( |l5| ) 
and others. Then we rediscover equation ( ^ ) for the 
coordinate. It implies that 



7 



2hti. 



(46) 



Equation for hr coincides with (Pq). Coercive field he is 
determined by equation: 



he = \/hti{2hQ - hti). 



(47) 



Eqn. (|41|) describes the shape of the hysteresis loop for 
ho > ht2 ■ It can be used also for the range of amplitudes 
hti < ho < ht2 on the lower branch of the hysteresis 
curve < h < ho- On the upper branch of this curve 
hr < h < ho the sign of the square root must be re- 
versed. 

The hysteresis loop area is determined by a following 
equation: 



A w 4:Mshr 



2Mshrjho 



X 2 



arcsm jr- 

hn 



v/1 - (hr/hor ^ho^\ 



(48) 
(49) 



This equation is valid for ho > ht2- If, on the contrary, 
ho < ht2, the area is: 



A = 2M,£^ [l-^\dh 



L 



2M. 



Jhr 



Z+{h) - Z_(/i) 



dh 



(50) 
(51) 



where Z± (h) denote the value of coordinate on the up- 
per (-I-) and lower (— ) branches of the hysteresis curve 
respectively. Performing the integration, we find: 

A = AM^K fl-l!^] + 7:Ahho^. (52) 
\ utL / ujL 

One can observe that for the case m ^ 1 the static thresh- 
old field hp does not enter physical results. Instead a new 
dimensionless parameter 



u = ujL/{"fho) = L/L^ 



(53) 



emerges. In the range of existence of the full hysteresis 
loop u' < 1/2. Therefore u' can not be large. However, it 



can be very small. It is worthwhile to consider separately 
the behavior of all values of interest in this asymptotic 
regime (large fields or small frequencies). For u' ^ 1 we 
find: 



he = hoVv/; 



A w AALK. 



(54) 



(55) 



(56) 



Near the threshold field hti the values of hr and A are 
given by following asymptotics: 



hr « 2^Jhtl{2ho - hti 
_ wMshti 



ATTM,^hti{2ho-hti) 



(57) 
(58) 



Near the second threshold field ht2 the reversal field hr 
is close to ho'- 



ho - 



{ho - ht2f 



2h 



t2 



(59) 



For m' ^ 1 we find scaling behavior of all hysteresis char- 
acteristics with universal critical exponents, independent 
on the pinning centers. Therefore, one can expect that 
the same scaling is valid for a clean ferromagnet with the 
relaxational dynamics. In particular, we have found that 



Ao^uj^l^hl'^ 



(60) 



In this section we assumed that the magnetization in a 
single domain reaches its saturation value. If it is not the 
case, the magnetization on the parts of hysteresis curve 
beyond the hysteresis loop ( "the whiskers" ) is not a con- 
stant. The exact shape of the hysteresis curve in this case 
depends on the single-domain equilibrium magnetization 
curve. 



IV. CREATION OF BUBBLES BY A MOVING 
DOMAIN WALL 

In this section we show that, at sufficiently high mag- 
netic field hp < ho < 77, the moving wall passes a series 
of defects which it cannot overcome. Instead it leaves 
closed domains of magnetization, opposite to the propa- 
gating magnetization, which serve as nuclei at the next 
half-cycle of the hysteresis process. 

So far we assumed that an almost planar domain wall 
is pinned by typical fluctuations of the random field ri{r). 
In a (spherical) region of the volume R^^d = D + 1 a. typi- 
cal fluctuation of J^^ ?7(r)d''r is of the order ±ri{l'^R'^y^'^ , 
the corresponding probability is of the order 1/2. In the 
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following we will consider the possibility of clusters of 
rare random field fluctuations which may serve as nu- 
clei of reversed spins when the field direction changes. A 
cluster of the size R consisting mainly minimal negative 
value —T] of the random field has the probability R) ~ 
exp {—{R/l)^ In 2). Such a cluster pins the domain wall 
for h < 7] even \i h > hp. In general, if the concentration 
of —rj sites in a cluster is c, its probability p{c, R) is of the 
order exp [-(i?/0''[cln c + (1 - c) In (1 - c) + In 2]] and 
the pinning force density is of the order (2c — 1)77. Since 
hp <C ?7, there is a field region hp < h < (2c— l)r] in which 
the domain wall can be pinned locally by these rare fluc- 
tuations of the random field. The true pinning condition 
is now given by (see also |2^ for similar considerations) 



r 

— + h-{2c- l)r] < 
it 



(61) 



where the first term denotes the Laplacian curvature 
force. Thus, the domain wall cannot overcome the pin- 
ning cluster if it has the size R bigger than Rmin, where 



(2c- l)r?- /i 



(62) 



Using ( |62| ) we get 

p{c,Rmin) ~ exp < - 



r 



(2c- l)r]l-hl 



5(c) (63) 



where 



g(c) = c In c + (1 - c) In (1 - c) + In 2 



(64) 



A closer inspection shows, that in c? > 2 dimensions 
p(c, i?inin) has its maximum at c — 1. The mean dis- 
tance between these strong pinning clusters is therefore 
of the order 



cluster 



I exp ■ 




(77 - ho)l 



Here we have replaced h by by the maximum field 
strength during one cycle ho- Since the interface can- 
not overcome the pinning cluster it will surround it and 
finally leave it behind as an island encircled by a domain 
wall (Fig. H). These islands serve as nuclei of the fa- 
vorite phase once the external field is reversed. Thus, in 
the expressions found in Sec. 



Ill 



the system size L has to 
be replaced by L duster as soon as i > L cluster- If 
approaches and finally exceeds ry, Lduster diverges and 
hence this type of magnetization reversal process disap- 
pears. Domain walls originate then either from surfaces 
or from a nucleation process which we consider in the 
next section. 



V. NUCLEATION CONTROLLED HYSTERESIS 

This mechanism works when the nucleation is the 
longest process. Let there exists a number of defects 
or impurities which lower the energy barrier for revers- 
ing magnetization locally. Still the reduced barrier A 
remains rather large so that the average reversal time 
Tr(0) — :/~^exp(A/T) {v is a microscopic frequency) 
is very large in comparison to the oscillation period: 
LJTr S> 1. However, the reversal time depends on mag- 
netic field: 



Tr{h) = Tr(0) BXp — 



Mgh 

T 



V ^ exp 



A - AL,h 



T 



(66) 



where Ma is the saturation magnetic moment in the ac- 
tivation volume. We assume also that lo/v <^ 1, but ac- 
cording to the accepted assumption exp(A/T) 3> v/oj. 
Then the probability of spin reversal in an individual nu- 
cleus is negligibly small if < - In (^27r^^ . It 

becomes reasonably large starting from the dynamical 
threshold for the amplitude /iq: 



^tl 



- Ti- In 27r- 



(67) 



Due to the sharpness of the exponential function, the re- 
versal time Tr{h) fastly becomes much smaller than the 
oscillation period 2ti/lo as soon as h{t) overcomes the 
threshold values hn- With the precision of a small pa- 
rameter jz/wexp (— A/T) the reversal proceeds at a fixed 
field h — hti which simultaneously plays the role of co- 
ercive field he and the reversal field h^ (see section |l|). 
With the same accuracy the reversal proceeds at a mo- 
ment of time tr = o;^^ arcsin (/iti/ft-o)- The hysteresis 
loop has a rectangular shape (see Fig||). With a little 
higher accuracy the transition proceeds at a value of the 
phase If — Lut determined by 



Mr,hr 



■ sm if 



A /2tti' 



(68) 



which results in a corrected value of the coercive field: 
hr — hr ^ hf^ — In ( arcsin -r^ ) (69) 



"*i ~ TT In arcsm — 
Ma \ ho 



As soon as the magnetization reverses in an individual 
nucleus, its walls expand fastly and, for a much shorter 
time than tr, the full magnetic moment reverses. The 
case in which the propagation of domain walls is the 
longest process has been considered earlier in Sections 
HI and IV. The HLA for a rectangular loop is simply: 



A 



4:Mshc 
Ms , 



T 27rzy T , 

— In —In 

A w A 



arcsm 



hti 
ho 



(70) 
(71) 
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Since Mahti ^ A the value 1 
For ho 3> hti the HLA reads: 



(r/A)log27ri//w < 1. 



^ = ^0 



1 



T , 2tti' 
— In 

A LU 



T. ho' 

^ 1 — 

A hti 



Note that at fixed lo and varying ho we find 



A^Ao{lu) 



Mahti hti 



(72) 



(73) 



and T/{Mahti) > T/A. This theoretical result can be 
compared with recent measurements by Suen and Ersk- 
ine [|l2| . They have measured the HLA for a thin film of 
-Fe/M^(110) in a wide range of frequencies u and ampli- 
tudes ho ■ They concluded that A ~ hgUj^ with small ex- 
ponents a ^ 0.2 and very small exponent /3 ^ 0.03 — 0.09. 
Both exponents depend on temperature. These conclu- 
sions are compatible with our result (|7|). Indeed, at such 
small values of exponents it is impossible to distinguish 
between the logarithmic dependence and the exponen- 
tial dependence with a small exponent. Our predictions 
a — T /Mahti and (3 — T/A agree with the experimen- 
tally discovered linear temperature dependence of the ex- 
ponent (3. It would be worthwhile to measure the fre- 
quency dependence of the dynamical threshold field hti- 
The fact of the existence of such a threshold has been 
observed in the same experiment. It should be noted 
that the value T /Mahti is not so small in the experiment 
~ 1/4. Therefore we expect that our theory is correct 
with the accuracy 25% only. On the other hand, the 
statement that there exists the power-like dependence of 
the HLA on ho is not proved convincingly by this exper- 
iment since the interval of variation for ho was too small 
(the amplitude changed about 10 times). 

So far we considered A and Ma as fixed values. This 
is correct if only one type of defects mediates the nucle- 
ation. More realistically these values are random. Then 
the defects with the minimal ratio A/Ma initiate the 
magnetization reversal process. 

According to Eqn.(|67|) hti becomes zero at very small 
frequencies w < 27: v exp (— A/T). Our theory is not valid 
for such a small frequency, but it reflects correctly the 
narrowing of the HL. On the other hand, the transition is 
smeared over the interval of magnetic field A/i ~ T /Ma- 
lt is small in comparison to h if hoMa 3> T. At a smaller 
amplitude ho the reversal is possible, but fluctuations 
of magnetization grow rapidly and become of the or- 
der of magnetization itself. In principle, it is possible 
that the nucleation time is of the same order of magni- 
tude that the time of domain wall propagation. However, 
since the nucleation time depends fastly on parameters, 
a small change of the regime makes one of the two times 
much larger than another. Thus, presumably the nucle- 
ation mediated hysteresis loop has the rectangular shape, 
whereas a curved shape indicates that the hysteresis is as- 
sociated with the domain wall propagation or with strong 
fluctuations of magnetization. 



A reason for the rounding of the HL in the case of the 
nucleation controlled hysteresis can be that the sample 
is splitted into a number of magnetically disconnected 
grains with different A and Ma in each grain. Then the 
shape of the hysteresis loop reflects the distribution func- 
tion for the ratio A/M^. 



VI. THE ADIABATIC DYNAMICS 

In the case of strong Ising anisotropy the dynamic 
equation for magnetization m in the continuum limit has 
the form: 



dm 



dm 



(74) 



where F is a free energy. In the adiabatic limit the total 
magnetization m{t) at the moment t can be represented 
as 



m{t) = mo(T, h{t)) + mi{t) 



(75) 



where mo{T, h{t)) is the equilibrium value of magnetiza- 
tion for a given momentarily value of the magnetic field 
h{t) and mi{t) is the deviation of magnetization from 
its equilibrium value. We discuss later restrictions im- 
posed by the adiabaticity condition. Substituting m{t) 
into Eq. (|74|) , one finds equation for mi: 



d^F 
dm^ 



mi, 



(76) 



where x = dm/ dh is the magnetic susceptibility. Eq.|7^ 
can be rewritten in a form: 



mi — h 

r 



(77) 



The hysteresis appears since the derivative h changes 
its sign when h passes its extremal values ±/io. For 
harmonically-oscillating field h{t) = /losinwt, its deriva- 
tive may be expressed in terms of the field itself: 



h = ±uJ^|Jll-hl- 



Plugging (^ into (||), we find: 



mi = T'^^hl-h^ 



(78) 



(79) 



Perturbation theory is valid if mi <C {mo)max — 
mo{T,ho), i.e. if ujx{ho)/T ^ 1. The coercive magnetic 
field he satisfies a following equation: 



mo(T, he) = \miiho,hc,T)\ ^^^h^^ hi (80) 
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Assuming he ^ ho, eqn. 
follows: 

mo{T,hc) = 



of) can be transformed as 



(81) 



ujx^ho 



The HLA A can be calculated by knowledge of mi: 



A 



4 / niidh 
Jq 



4lu 



hi 



h^dh. 



(82) 



Near Curie point the relaxation time Tc oc T^^x grows as 
Tc ^ T^e^'^^ where e = (T — Tc)/Tc, v is the correlation 
length critical exponent {v — 1 for 2D Ising model) and 
z is the dynamical critical exponent. The best numeri- 
cal value for z is z ~ 2.17 (for a modern review see ||2^ ] 
and references therein). If e = and h ^ 0, the relax- 
ation time is proportional to h^^^^^. For adiabaticity it is 
necessary that uj <$:. ujc — 'max{TQ^e^,TQ^{hQ/hex)^''/^^), 
where h^^. is the saturation (exchange) field. 

The amplitude /ig should be considered as small (the 
case i) or large (the case ii) depending on its ratio to 
a characteristic field h*{e) = he^^^^^^- For a small am- 
plitude ho <C h*{e) one can employ linear approxima- 
tion for mo{T, h) = x(e)/i to solve eqn.(|l]). Then he = 
/iowx(e)/r(e). Next we use the relation: xi^)/^ tqe^'' 
to find: 



ujhQTQe ^ = ujhoToe 



-2.17 



(83) 



For a large amplitude ho S> h*[e) a dimcnsionless ratio 
K — ti;x(e)/r(e) matters. If it is much less than the large 
value ho/h*, then the coercive force can be calculated as 
he = Kho, exactly as in eqn. (^). Otherwise h^ 3> h* . 
Then an approximate solution of eqn. (^l|) reads: 



he ~ hex ^To 



LUTo-. 



-0.46 



hex J \ h 

Corresponding results for the HLA are as follows: 



(84) 



i) Acc ujhQ- 



uh^Toe 



uhlroe-^-^^. (85) 



ii) A cx Luho 



h*{e)x\e) 



LohohexTo^ 



(86) 



Though we can not calculate explicitly the hysteresis 
loop area in the opposite, anti-adiabatic regime LUe ^ 
Lu <C J/h, the dimensionality arguments lead to a nec- 
essary estimate. Indeed, according to general equation 
(p2|), its scaling dimensionality is 



A(x xhl, 



(87) 



where we have used the facts that uj and T/x have 
the same dimensionality and that ho <C h*{Lo) — 
hex{huj/jf^'^. Thus, 



A oc x(w)/io oc UJ ^^^h^ = u) 



-0.8l2 



Typical orders of magnitude are: Tq ^ ~ J/h ^ 4 • 
10"Fz;e - l{)-^]UJe{e) - 2MHz;h*{e) - 30e;Te - 
300K; T~Te^ 0.3K. 

Note that all calculations in this section have been 
made for T > Te, so that we could ignore the problem 
of nucleation and domain wall motion. Below the Curie 
temperature the problem of spinodal decomposition ap- 
pears, when the field reaches the spinodal. Then, again, 
either the nucleation or the domain wall motion prevails 
in the magnetization reversal. It docs not happen if ei- 
ther ho S> h*{e) or w ^ T~^. Then one can use previous 
results. However, if these conditions are not met, then 
either nucleation or the domain wall propagation pre- 
vails in the magnetization reversal process. To establish, 
what is more important, let us consider corresponding 
time scales. The nucleation time for a pure sample is: 



Tr = Te{e) ■ exp 



h*{e) 
h 



(89) 



Though this time scale grows at approaching to the Curie 
temperature, it can be sufficiently small. Its minimiza- 
tion over e at fixed value of h gives: 



Tr OC To 



^/(/3+7) 



To 



''CX 



1.157 



(90) 



For the exchange field about 10^ Oe and external field 
about lOOOe the nucleation time becomes of the order of 
10~^sec. Thus, the nucleation proceeds very fastly and 
the process of domain wall propagation prevails in this 
range of fields and frequencies. This conclusion explains 
the results of our numerical calculations for temperatures 
close to the Curie point (see section VII). 



VII. SIMULATION OF THE DOMAIN GROWTH 

In this Section the process of magnetization reversal 
in a few monolayers thick magnetic film is modeled by 
analogous process in the two-dimensional Ising model 
with the Glauber dynamics. In simulation we neglect 
the dipole-dipole interaction and demagnetization effects. 
The Hamiltonian of the Ising system in a time-dependent 
field h{t) has a form: 



(91) 



where summation over r runs over the lattice sites and 
a labels nearest neighbors. In a perfect ferromagnetic 
film Ja(r) = Jo > independently on r. In the case of 
quenched disorder we assume that the exchange integral 
is ferromagnetic (Ja(i") = Jo > 0) with the probability 
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1—1? and antiferromagnetic (Ja(r) = — Jo < 0) with the 
probability ( d <^ 1). No correlation in the location of 
random bonds is assumed. In Sec.|| it is shown that this 
type of disorder is of general importance for the problem 
of the domain wall moving in the disordered media. Weak 
random bond disorder does not destroy the long-range 
order in the 2D Ising model. The magnetic field h{t) is 
supposed to oscillate harmonically as h(t) = /iq sin (wt), 
unless a different assumption is not specially formulated. 
The process of magnetization reversal can be divided into 
two stages: nucleation of domains with opposite mag- 
netization and growth of these domains. Which process 
dominates depends on system parameters and its history. 
In this Section we consider a limit of the nucleation time 
much smaller than the growth time. Then the magneti- 
zation reversal process is dominated by the domain wall 
propagation discussed in Secjl| and [II. 



We employed the Monte-Carlo simulation with the 
Glauber dynamics (see e.g. |Q) to check Eq.^for perfect 
systems and systems with disorder. We modelled the dis- 
order by a small concentration of randomly distributed 
antiferromagnetic bonds. This disorder weakly influences 
the phase diagram and results only in a small shift of Tc, 
where Tc is an Ising transition temperature of the model 
Eq.^ with h = 0. We have found the linear dependence 
of domain wall velocity on applied magnetic field at all 
relevant temperatures with and without disorder. Ap- 
parently, in these simulations we dealt with the limit of 
weak pinning, so that hp <^ Hq. However, at low temper- 
atures (T < 0.1 • Tc) even a small disorder, which almost 
does not change the thermodynamics, modifies the do- 
main wall dynamics drastically. The results related to 
the domain wall motion and hysteresis phenomena in a 
weakly disordered system at low temperature will be dis- 
cussed elsewhere. 

The domain growth was studied in two geometries: 
stripes and circles. To distinguish the domain growth 
from domain nucleation we examined either a specially 
prepared defect for fast nucleation, or a prepared nucleus 
with the opposite direction of magnetization possessing 
the shape of the circle or stripe. In both cases the magne- 
tization changed in accordance to the model of a straight 
DW, i.e. similar to Eq.^ 

Simulated hysteresis loops for different values of mag- 
netic field amplitudes {ho = 0.025J, ho = 0.05J, and 
ho = 0.125J) and constant period of 8 • 10"^ Monte-Carlo 
steps are shown in Fig. Hysteresis loops for differ- 
ent periodes (4-103, 8 • 10^ and 256 • 10^ Monte-Carlo 
steps) and constant ho — 0.05 J are shown in Fig.|[ The 
observed hysteresis loop behavior obey the general clas- 
sification scheme developed in Sees. ||,|l| (compare with 
Fig. |2|). In particular, the asymmetric loops in Figs. ^,|^ 
correspond to the loop in Figj^a, the loops with whiskers 
in Figs. |4||^ corresponds to the loop in Fig.^, and the 
symmetric loops without whiskers in Figs. ^,|^ corre- 
sponds to the loop in Fig.||c;. The concentration of ran- 



dom bonds in the system was « 0.8%. 

Fig.^ show the time sequence of inflating domains for 
the system with T = O.ST^ and ho = 0.25J. The con- 
centration of random bonds in this case was « 3%. The 
nucleation first proceeds at three regions, then the corre- 
sponding domains inflate rapidly and two of them merge. 
The time development of the inflation process is shown in 
Fig.^. The total magnetization for each snapshot is indi- 
cated in the Fig. ||. The analysis show that the inflation 
of a domain with reversed magnetization is described by 
simple law Eq.p|. 

Let Lat be the mean distance between nucleation cen- 



ters. It plays the same role as the system size L in Sec.|H. 
The ratio Lj^ / = ujv is an analog of the dimensionless 
parameter u' introduced by Eq.M. In the limit L^r <C 



(compare Sec. [II) the domain growth will be fast enough 
to provide the full hysteresis cycle with hr <^ ho (see 
Fig.^ ). We evaluated the scaling behavior of he and 
^ as a function of w and ho in the limit Lat <C us- 
ing Eq.^. To deal with the domain wall propagation 
dominated hysteresis, we considered temperature close 
to Curie point. We argued earlier (see Sec. VI) that in 
this case the nucleation proceeds lastly enough. Fig.|^ 
shows domain distribution in the applied reversed mag- 
netic field at r = 0.95Tc. Multiple nuclei which are seen 
on this picture serve as evidence that the nucleation is 
really fast enough. 

At h = he the typical domain size is of the order of 
the mean distance between domains. Since Ln <^ we 
conclude that he <C ho, and Ln ~ {"f / houj){hc)'^ . Assum- 
ing that the average distance between nucleation centers 
does not depend significantly on w, we immediately ob- 
tain that the coercive field he scales as he ^/uj (com- 
pare with Eq.^4|). If, in addition, the average distance 
between nucleation centers does not depend significantly 
on ho, we have he oc ^/haut. If we assume that the mag- 
netic field is strong enough and saturation magnetization 
does not change significantly by changing ho, the HLA 
area scales in the same way a.s he, i.e. A oc \/hoUj. This 
is a simplified qualitative version of arguments leading to 
Eq.|^ of Sec. 



[II 



Fig.|8| shows the dependence of the hysteresis loop area 
on w in the double logarithmic coordinates at T = Tc- 
The ^Juj scaling behavior is valid over three decades of 
the oj variation. However the the range of magnetic field 
values available is only one decade. For this reason we 
have only checked that magnetic field dependence is con- 
sistent with A cx \/ho dependence. 

To check the proposed mechanism of magnetization re- 
versal we have also simulated this process for the mag- 
netic field h = hosiv?{ujt). In this case both he and A 
scales in accordance with the analytical results as w'^/'*. 

Though we put T = Tc, the system is far away from 
the critical point due to relatively large magnetic field. 
Indeed, one can estimate e defined in Sec.^^ from the re- 
lation e « {ho/he^f^^^ (see Sec.0). In our simulations 
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ho/h^^ = 0.05-0.005 which give e « 0.2-0.06. It means 
that the deviation from the critical point is large enough. 
With these values of fields our system is also far away 



from the strong fluctuation regime described in Sec.V] 



A necessary condition for this regime is that mo ^ 1. 
The magnetic moment mo scales as mo c>c (h/hexY^'^^ ■ 
Therefore the field ft-o ~ O.OOS/iea; is rather strong it gives 
mo ~ 0.7. 

With the temperature decreasing, away from critical 
region the HLA deviates strongly from the power-law be- 
havior as it is follows from Fig.|[ We expect that the HL 
changes its shape and characteristics from those domi- 
nated by domain wall propagation to those dominated by 
nucleation. The crossover behavior can not be described 
by a single power law. 



VIII. CONCLUSIONS. 

We studied main mechanisms of the hysteresis in a 
ferromagnet: the driven DW motion, nucleation and re- 
tardation in non- linear magnetization dynamics. Each of 
them may be dominant at proper conditions. A process 
is dominant if it provides maximum value for the coercive 
field he- 

On the background of general theory of the interface 
motion in a random medium we studied the hysteresis 
process controlled by the DW motion. We introduced 
two dynamic threshold fields hti and ht2 corresponding 
to the occurrence of the full magnetization reversal and to 
the occurrence of the single-valued parts on the hystere- 
sis curve respectively. These dynamical threshold fields 
are larger than the static threshold field hp which is re- 
quired to start a motion of the DW. We established that 
hti and ht2 measured in units hp are functions of one 
dimensionless parameter u — ujL/^hp (see notations in 
sect. I, 1^ HI). The coercive field he and the reversal 
field hr, measured in the same units, can be expressed in 
terms of one function of two dimensionless parameters: 



hr. = hr,F 



LUL hr. 



2-fhp ' hp 



Experimental observation of this type of scaling would be 
the best indirect evidence of the DW motion controlled 
hysteresis. The direct observation of the DW motion by 
SMOKE and STM methods in principle is possible |2|] 
and is much appreciate. At large fields /lo ^ hp the de- 
fects are inessential. Therefore, the dependence on hp 
must vanish from scaling laws (|9^). It happens indeed, 
and both fields he and hr are expressed in terms of one 
dimensionless parameter: 



where Fix) = 2yjx{l — x). Note that in this limit the 
value hp vanishes from equations for the dynamic thresh- 
old fields as well: hti — ^-^/7, ht2 — 2/iti. We pre- 
sented also corresponding equations for the HLA to which 
the most experimental efforts were concentrated. How- 
ever, we would like to emphasize that the HLA is not 
the only measurable characteristics of the HL and even 
not the most informative its characteristics: the fields 
hti, ht2, he, hr as well as the shape of the hysteresis curve 
are not less interesting. The functional dependence of 
the fields hti and ht2 on ho at ho not much overcoming 
the static threshold field hp allows to restore the basic 
function f{{h — hp) /hp) of theory of the DW motion in 
a defect medium. 

For sufficiently large ho the length L in equations ( |9^ , 
|9^ ) and others is the size of the system. Thus, the size 
effect is observable in the HL characteristics. However, 
at fields, larger than hp, but still smaller than the maxi- 
mum defect strength, the moving DW produces bubbles 
playing the role of ready nuclei for magnetization reversal 
at the next half-cycle of the hysteresis loop. Thus, the 
length L in this case is the average distance between the 
bubbles. Its strong dependence on the the field ampli- 
tude ho makes the scaling laws ( |92| , p3|) less transparent, 
but they become size- independent. The scaling becomes 
simple again if the dominant defects are of topographic 
origin. The density of topographic defects is independent 
on the amplitude of magnetic field. Then L is the average 
distance between such defects. 

In the case, when the driven DW are almost free 
{ho ^ hp) and the HL is narrow {ho ^ hti), the HLA 
was found to be proportional to w^^^/io^^. This conclu- 
sion is supported by our numerical MC simulation. 

In the nucleation controlled process almost rectangular 
HL shape is expected, unless the sample is divided into 
a multitude of magnetically disconnected grains. Note 
that the shape of almost vertical parts of the HL in this 
case is determined by the DW motion and, therefore, is 
universal. In the case of nucleation controlled hysteresis 
the HLA must grow logarithmically with the frequency 
and amplitude of magnetic field. However, the coefficient 
at the logarithm of frequency is smaller than that at the 
logarithm of the amplitude. Slow nonlinear critical dy- 
namics may produce the HL with various scaling limits 



as described in sect.V] 



he = hoF 



hr = hoF — — 



(93) 



The comparison with the experiment is still rather 
poor since there is no systematic study of the HL char- 
acteristics as functions of dimensionless parameters. The 
experimental efforts were focused on the HLA, but the 
precision of the measurements is not high enough for re- 
liable determination of scaling. Sue and Erskine re- 
ported the observation of linear dependence of the HLA 
on the logarithm of frequency ]30[ | with the coefficient 
proportional to temperature. Both these facts agree with 
our theory. However, they also claimed the power depen- 
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dence of the HLA on the amplitude ho with the exponent 
about 1/4. This dependence contradicts to our theory. 
The reason for discrepancy is not yet clear. The interval 
of variation of Hq is not large enough to establish the ex- 
ponential dependence reliably, but it definitely deviates 
from logarithm. On the other hand, our theory gives 
logarithmic dependence on ho only if ho ^ /iti- Further 
studies both theoretical and experimental are necessary 
to clariiy the situation. 

Our numerical simulations show visibly the formation 
and propagation of domain walls as the dominant pro- 
cess at high temperature. The domain walls look rather 
rough. It can happen as a consequence of strong criti- 
cal fluctuations which are not taken into account in our 
DW motion theory, or it may follow from the fact that 
the dimensionality 2 is the marginal for the development 
of roughness due to defects (the roughness exponent <^ 
equal to 1 for d — 2). Nevertheless, the scaling low 
A cx ' /iq for h ^ hp is confirmed in numerical sim- 
ulations at high temperatures. At lower temperature the 
nucleation becomes dominant. 

An interesting possibility to verify our predictions for 
the domain growth controlled hysteresis is a study of a 
hysteresis of adsorption isotherms in a close to equilib- 
rium conditions. The process of close to equilibrium ad- 
sorption can be studied experimentally for noble gases 



this work has been completed. He also thanks Prof. J. 
Zittartz for the hospitality extended to him at Cologne. 



(see e.g. |2^]). The quasiec 
films has been studied in 



uilibrium desorption of helium 
2^ . In the case of equilibrium 
adsorption the chemical potential plays the role of mag- 
netic field and the coverage - plays the role of magnetiza- 
tion. The domain wall width for adsorbed systems is only 
one lattice period wide. For this reason one can expect 
rather small nucleation barriers for equilibrium adsorp- 
tion and more broad range of existence for domain growth 
controlled hysteresis than in magnetic films. More de- 
tailed analysis is given in ||2^. Though the adsorption 
isotherms have been measured starting from Langmuir, 
no hysteresis loop measurements has been performed un- 
til recently. The first such an experiment by H. Pfniir 
and K. Budde 129] gave values of P close to 1/2. 
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FIG. 1. Schematic picture of the ratcheted potential in the 
random bond (a) and random field (b) case. 



FIG. 2. Schematic pictures of hysteresis loops. 

(a) Incomplete HL for ho < hti. 

(b) Symmetric HL for ho — hti. 

(c) The HL for hti < ho < ht2. 

(d) The HL for ho > ht2. 

The values hp, he, hr, ho are marked in all figures. 



FIG. 3. Formation of a bubble by the moving domain wall. 



FIG. 5. Simulated hysteresis loops for different values of u) 
and constant ho- The observed hysteresis loop behavior obey 
the general classification scheme developed in Secs.|^,|l^ and 
shown in Fig. S (see explanation in text). 



FIG. 6. Infiation of domains with opposite magnetization 
in an Ising model in reversed magnetic field with ho — 0.25J 
at T = 0.5Tc for the system of the size 192 x 192 with pe- 
riodic boundary conditions and with random frozen defects 
(antiferromagnetic bonds) with concentration 0.03. Each dot 
corresponds to the bond with opposite spin orientation. 

FIG. 7. Domain walls distribution in an Ising model in 
reversed magnetic field with ho = 0.05J at T = 0.95rc for 
the system of the size 432 x 432 with periodic boundary con- 
ditions. Each dot corresponds to the bond with opposite spin 
orientation. 

FIG. 8. Hysteresis loop area as a function of frequency in a 
double logarithmic coordinates for the magnetic field with am- 
plitude ho = 0.05 J at T = Tc (lower line) and at T = 0.95rc. 



FIG. 4. Simulated hysteresis loops for different values of 
ho and constant u. The observed hysteresis loop behavior 
obey the general classification scheme developed in Secs.|l[ 
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and shown in Fig. ^ (see explanation in text). 
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